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Abstract. The presence of pronounced electronic correlations in a one-dimensional 
systems strongly enhances Coulomb coupling and is expected to result in distinctive 
features in the Coulomb drag between them that are absent in the drag between two- 
dimensional systems. In this article, we review recent Fermi and Luttinger liquid 
theories of Coulomb drag between ballistic one-dimensional electron systems, alias 
quantum wires, in the absence of inter-wire tunneling, to focus on these features and 
give a brief summary of the experimental work reported so far on one-dimensional drag. 
Both the Fermi liquid (FL) and the Luttinger liquid (LL) theory predict a maximum 
of the drag resistance Rd when the one-dimensional subbands of the two quantum 
wires are aligned and the Fermi wave vector kp is small, and also an exponential decay 
of Rd with increasing inter- wire separation, both features confirmed by experimental 
observations. A crucial difference between the two theoretical models emerges in the 
temperature dependence of the drag effect. Whereas the FL theory predicts a linear 
temperature dependence, the LL theory promises a rich and varied dependence on 
temperature depending on the relative magnitudes of the energy and length scales of 
the systems. At very low temperatures, the drag resistance may diverge due to the 
formation of locked charge density waves. At higher temperatures, it should show 
a power-law dependence on temperature, Rd oc T x 1 experimentally confirmed in a 
narrow temperature range, where x is determined by the Luttinger liquid parameters. 
The spin degree of freedom plays an important role in the LL theory in predicting the 
features of the drag effect and is crucial for the interpretation of experimental results. 
Substantial experimental and theoretical work remains to be done for a comprehensive 
understanding of one-dimensional Coulomb drag. 
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1. Introduction 

Moving charge carriers in a conductor exert a Coulomb force on the charge carriers in a 
nearby conductor and induce a drag current in the latter via momentum transfer. This 
phenomenon, known as Coulomb drag, was predicted by Pogrebinskii in his pioneering 
paper |jj in which he argued that in a structure of two semiconductor layers separated by 
an insulating layer, there would be a drag of carriers in layer 1 ("drag layer"), resulting 
in a drag current Id, due to the direct Coulomb interaction with the carriers in layer 2 
("drive layer"), where an electric current I flows. If no current is allowed to flow in the 
drag layer, the charge carriers will accumulate at one end inducing a charge imbalance 
across the layer. This charge will continue to accumulate until the force of the resulting 
electric field balances the drag force. In the stationary state there will be an induced 
or 'drag' voltage Vd in the drag layer. When the carriers in both layers are of the same 
type (electrons or holes), the drag voltage has a sign opposite to the voltage drop in 
the drive layer. Figure [I] gives a schematic view of Coulomb drag between two parallel 
quantum wires. The quantity usually measured in experiments is the drag voltage Vd- 
The drag resistance Rd is defined as Rd = -Vd /I. 



Figure 1. Schematic view of Coulomb drag between parallel quantum wires. 

Coulomb drag between two-dimensional (2D) electron systems has been extensively 
studied [|| both experimentally and theoretically. The basic physics involved in the 
description of the drag in two dimensions is now well understood on the basis of 
Fermi liquid (FL) theory of interacting fermions. The FL theory is well established 
in three dimensions, holds marginally for many two-dimensional systems, but generally 
fails in one-dimension. The theory is based on Landau's conjecture that the low-lying 
excitations of interacting fermion systems can be connected continuously to those of the 
non-interacting Fermi gas - there is a smooth mapping between the quasiparticles of 
the interacting and of the non-interacting system 0. 

Coulomb drag between one- dimensional (ID) electron systems has been the focus 
of considerable interest in recent years because our understanding of the quantum 
properties in interacting ID systems is unsatisfactory. Experimental work on the subject 
remains quite limited [|], [| ||, |7j , a fair number of theoretical papers have been published 
H |9], [T(J [11], [12], [L3|, [14], [L5|, [16], [L7|, [18j. The primary reason for this theoretical interest 



is that Coulomb drag is one of the most effective ways to study electron-electron (e-e) 
interaction. 

It is now theoretically established that in an interacting ID electron gas of 
infinite length the e-e-interaction completely modifies the ground state of the system. 
The elementary excitations can not be treated as non-interacting quasiparticles of 
a conventional Fermi liquid, but instead acquire a bosonic nature. An adequate 
theoretical description of these interacting ID systems can be done in terms of the 
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so-called Luttinger-liquid (LL) |19| (for a recent review see pOfl), complementary to 
the FL description in higher dimensions. In a real ID system of finite length at finite 
temperatures, the extent of influence of the e-e interaction will depend on the system 
parameters. 

Experimental efforts to observe manifestation of Luttinger liquid behavior, however, 
have been quite limited PIL f22|, ^ . Part of the problem is that the (e-e) interaction 
has little influence on the conductance of a single wire, since the current is proportional 
to the total electron quasi-momentum, which is conserved in electron-electron collisions. 
To look for experimental evidence of the LL state it would therefore help to explore 
a new experimental tool based on new devices and physical phenomena. Coulomb 
drag between ID electron systems in a dual-wire configuration and opens up a new 
opportunity and avenue for experimentally probing the LL state in a ID electron system. 

The purpose of this review paper is to present the current status of the theory of 
Coulomb drag between ID electron systems for electron transport in the ballistic regime, 
and to report on experimental measurements of the ID drag effect. Ballistic transport 
takes place when the quantum wire dimensions are smaller than both the elastic and 
the inelastic scattering lengths. Electron transport is strictly one dimensional when 
only the lowest ID subband of the wire is occupied and transport takes place in the 
fundamental mode. The ballistic regime is well suited for Coulomb drag study, since 
in this regime other scattering processes, such as impurity and phonon scattering, are 
either insignificant or totally absent. 

Further theoretical and experimental investigation of the ID Coulomb drag 
effect can enhance our general understanding of the properties of systems of low 
dimensionality. This broad class of systems is currently a very active area of research. 
In addition to its fundamental interest, a comprehensive understanding of Coulomb 
interaction between quantum wires is expected to play a significant role in the design 
of nanodevices, such as single-electron transistors (SETs) p5j and quantum cellular 
automata (QCA) f26j, which are comprised of quantum dots and quantum wires in 
close proximity. 

The paper is organized as follows. Section describes the theory of ID Coulomb 
drag based on the Fermi liquid approach. Section || reviews the Luttinger liquid 
description of this effect. In Section [| is presented a summary of the experimental 
work reported so far and a comprehensive analysis of the experimental results using 
both the FL and the LL descriptions of ID Coulomb drag. Finally, Section [5] gives some 
guidelines for future work on the subject. 



2. Fermi liquid approach 

In this section the ID Coulomb drag is analyzed within the Fermi-liquid concept. We will 
follow Refs. H, [T8| and use the physical picture developed by Landauer P8[ , Imry P§ 



and Buttiker f3(J. We assume that each quantum wire is connected to ideal electronic 
reservoirs attached to its ends. The relaxation processes in the reservoirs are considered 
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to be so fast that each of them is in thermal equilibrium. 

The e-e interaction within a single quantum wire does not result in a current 
variation because of the quasimomentum conservation in the e-e collisions. However, if 
two such wires, 1 and 2 are near one another and are parallel, the Coulomb interaction 
of electrons belonging to different wires can transfer momenta between the wires, which 
eventually gives rise to a drag effect. 

The drag force due to the ballistic current in wire 2 creates a sort of permanent 
acceleration on the electrons of wire 1. As wire 1 has a finite length L a steady drag 
current J is established. 

Within the Fermi liquid approach we restrict ourselves to direct electron-electron 
collisions mediated by the Coulomb interaction. Let us analyze the conservation laws 
for such collisions of electrons belonging to two different wires, 1 and 2, each of them 
being parallel to the x-axis. We have 

e W+ e {2) -e (1) +e (2) (I) 

t nk * t n'k' ~ H,k+q ' t l',k'-q \ L I 

Here hk is the x-component of the of the electron quasimomentum. In this and the next 
sections we will put fi = 1 and k& — 1 (where fee is the Boltzmann constant). These 
quantities will be restored only in some final (or most important) formulas. Now, 

e™=eW(0)+k 2 /2m, (2) 

m being the effective mass while n being the transverse quantization subband (channel) 
index, with primed quantities corresponding to wire 2 throughout. The solution of Eq. 
(H) can be written as 



q = -{k - k')/2 ± y/(k - k'Y/A + mbe (3) 

with 6e = e£\0) + e { *)(0) - e, (1) (0) - 6^(0). 

We assume that the electrons of the quantum wires are degenerate and the 
temperature is low compared to the electron Fermi energy. For the electron-electron 
collision to be possible, the absolute values of the four quantities, namely, e^, e^) k ,, ep fe +g> 

(2) 

and e\, k ,_ should be within the stripes k^T near the corresponding Fermi levels. This 
means that within the accuracy mT/kp the following relations should be valid 

u _ u _ £>') 



\k + q\=k%\ \k'-q\ = k [ p > . (4) 

Here k^ denotes the Fermi quasimomemtum for band n. In general it is impossible 
by variation of a single quantity, i.e. the transferred quasimomentum q, to satisfy both 
relations of Eq. (f|) (provided, of course, that the distances between the channel bottoms 
are much bigger than T.) 

In other words, one cannot in general satisfy Eq. (|l|) for a finite 5e. Therefore for 
a general case one should have n = I, n! = I'. If both wires are identical equations 
n = I', n' = I are also possible. In both cases Se = 0. We will assume the wires to be 
different. Then 

+ 41 - - 4 2 L 9 ) = (m/\q\)8(k -k' + q). (5) 
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This means that the quasimomentum transferred during a collision is q = k' — k, i.e. 
the electrons swap their quasimomenta as a result of collision. 

Assuming that the drag current in wire 1 is much smaller than the ballistic current 
in wire 2, we calculate the drag current by solving the Boltzmann equation for wire 1 
(otherwise we should have solved a system of coupled equations for both wires). We 
assume the wires to be different though having the same lengths L and consider the 
interaction processes when electrons in the two quantum wires after scattering remain 
within the initial subbands e£2 = £n\o) + k 2 /2m and e%) k = £^(0) + k 2 /2m, n being 
the subband's number. The Boltzmann equation for the electrons occupying the nth 
subband is 

v k ^=l^{F^,F^} (6) 

where i^ 1 ' 2 ) are the electron distribution functions in wires 1 and 2 respectively, and / 
is the collision integral. We assume that the only type of collisions that is essential is 
the interwire e-e collisions described by the term 

x (i2)| F (i) ;F (2)-j. = 2 Jdk^Jdq 

k + q,n; 2, k' — q,n'<— 1, k, n; 2, k', n')V (7) 



where 



P (l)p(2) (, _ \ (-, _ p (2) 

r nk r n'k'\ L 1 nk+q ) y L r n'k>- 



^nk+q^n'k'-q ^nk^j ^i'fc') ■ (^) 

2 is the spin factor; the scattering probabilities are assumed to be spin-independent. 
If the e-e collisions can be treated within the perturbation theory then the scattering 
probability is given by 

w(l, k + q,n,2, k' — q,n' <— 1, k, n; 2, k', n') = 
2tt| (1, k + q,n; 2, k' - q, n'\V\l, k, n; 2, k\ n')\ 2 x 

Af^ 1 ) j_ ,( 2 ) _ J 1 ) _ ,( 2 ) ^ (a) 

°\ t nk ^~ t n'k' t n,k+q t n',k'-qJ- K*) 

The matrix element of electron-electron interaction can be transformed to 

(1, k + q,n; 2, k' - q, n'\V\l, k, n; 2, k', n') 

= \J d 2 r ± J rf 2 r^;(r ± )0;,(r' ± ) x 

V q (r ± -r' ± )<i> n (r ± )<i> n ,(T' ± ) (10) 
where V q — J dxV(x, r±) exp(— iqx), = (y,z). We have 

J dx J dx' V(r - r')e iq{x - x,) = 2e 2 LK (\q\ \Ar±\) (11) 
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where Ar_|_ = r_|_ — r'j_ and K is a modified Bessel function defined in Ref. [331. Now 



It means that the e-e interaction goes down exponentially provided |?||rj_ — r'jj/ft 3> 1. 

To calculate the current in wire 1 we iterate the Boltzmann equation (^) in 
the collision term T^ 12 \ The first iteration gives for the nonequilibrium part of the 
distribution function AFn 



np 



^) -_( z± £] lj(l2) {j p(D ? p(2)y (13) 
\ 2 J v n 

for k > (fc < 0) respectively. One gets for the drag current 



id 



- 2eL E/°°§ J(12) ^ (1) ^ (2) }- 



(14) 



We assume in the spirit of the Landauer-Buttiker-Imry approach the driving wire 
connected to the reservoirs which we call "left" (+) and "right" (— ), each of these 
being in independent equilibrium. Let the x-component of the quasimomentum of an 
electron in wire 2 before scattering be k', after scattering by an electron of wire 1 be 
k' — q. Let k' > while k' — q < 0. Then the first distribution function in wire 2 
* s = f^'k ~ A*^) where / is the equilibrium Fermi function. The second one is 
^i\l+q = f{ e v,k+q~ Z^ - ^) where = fi±eV/2. At eV = the wires are in equilibrium. 
We denote the corresponding equilibrium chemical potential as \i. 

Let us denote by A{F} the expression one gets after substitution of the equilibrium 
distribution functions given above into the collision term. For k' > (k' < 0) and 
k' — q < (k' — q > 0) where k! is the electron quasimomentum in wire 2 before the 
scattering we obtain 

A{F (D )F (2) } = ±2sinh 

x^-f(^l-n)][l-f(e%,-^)} 

^f(^ k+q -n)f(e%,_ q -^). (15) 

We begin with a discussion of the Ohmic case eV/T <C 1. Accordingly, we replace 
sinh(eVy2T) by its argument and all the chemical potentials in Eq. (|T5| ) by the same 
value /i. The initial and final states of the colliding electrons should be within T of 
the Fermi levels ||. This means that only the terms with e« (0) = £^(0) where the 
equality is satisfied with the indicated accuracy give the principal contribution to the 
current. (The importance of equal channel velocities was also pointed out in Ref. ||). 
The contribution of each such pair of levels to the current is 

_ e 5 m 3 Lk B TeV 1 , [e£ } (0) - 



X 



2^ kn ^nny^n, 

eS\0)-e%>{0)' " 



sinh 



2k B T 



(16) 
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where 



9nn'{q) 



d 2 r ± / d 2 r' ± \<t> n (r±)\ 2 \Mr±)\ 2 K (q\Av ± \) 



(17) 



K is the dielectric susceptibility of the lattice, k n = ^mj/i — £ n (0)]. This equation 
has been also obtained using linear response theory. In the case considered, wire 2 is a 
part of a usual structure for measuring ballistic conductance, i.e. it joins two classical 
reservoirs, each of them being in independent equilibrium. The driving current is 

I = M-V, (18) 

7T 

M being the number of active channels (i.e. the subbands whose bottoms are below the 
Fermi level). So far a simplifying assumption has been used: the chemical potential fi 
in wire 1 and the average chemical potential in wire 2 are equal. In the general case 
they can have different values fi^ and respectively. Then one still gets Eq. ([T6|) 
with the replacement 



One can measure either the current or the voltage that builds up in wire 1. The 
ratio of the drag current to the ballistic driving current for in\o) = e^, (0) is given by 

in TVrV'K J\l 



where 



Vnn> = — — 3 W (tf 5 + fcW) . (20) 

( k n -\- kn 



Here kn )2 ^ = y 2m[/x^ 1 ' 2 ) — En '^(O)]. In this approximation, k n 1,2 ' > = k^,' 2 \ In an 
experiment one usually measures the drag resistance Rn = —Vd/I = I^Gd/I, where 
Gd is the ballistic resistance of the drag wire and depends on the number of occupied 
subbands. 

The ID subband structure of the wires can be modified by changing the effective 
wire widths by applying appropriate gate voltages (Fig. §). The variation of gate 
voltage may affect the positions of the levels of transverse quantization in the two wires 
in a different way. In the course of such a variation a coincidence of a pair of such 
levels in the two wires may be reached. The estimate (|20D is not very sensitive to the 
form of confining potential and electron densities. In Fig. |3] the ratio I /Id is plotted 
(for yS^ = /i*- 2 -*) as a function of the ratio of effective wire widths. This plot exhibits 
striking oscillations with large peak-to-valley ratios. The peaks occur when channel 
velocities in two interacting wires are equal which happens whenever any two current- 
carrying channels line up. This sort of coupling is particularly strong when such channel 
velocities are quite small. 
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The condition ei^O) = e£)(0) gives the main maxima of the drag current, especially 
for the lowest levels of transverse quantization. Some subsidiary maxima can also be 



observed, particularly in external magnetic field — see Ref. [14 . 

So far we have considered the peaks of the drag current under conditions where 
the Fermi level is well above the coinciding bands in wires 1 and 2. Now we would like 
to say a few words about a special case that may be particularly important regarding 
the experiment described below. This is the case where the bottoms of two subbands 
not only coincide but also just touch the Fermi level. Then the conduction electrons 
obey the so-called intermediate statistics. It means that their equilibrium distribution 
functions are 

1 

exp(k' 2 /2mT) + 1 

The drag current is proportional to the e-e scattering probability averaged over this 
distribution function. The scattering probability itself is determined by the quantum 
mechanics, i.e. it depends on the form of electron wave function in the quantum 
wire, in other words, on the exact form of the confining potential. Investigation of 
the temperature dependence of the drag in this special case is one of the problems of 
the theory to be solved in future. 

Now we turn to the case of non-Ohmic transport in the drive wire, i.e. to the case 
where eV ^> k^T [181. The situation for /i« = ^ is illustrated in Fig. ||. The upper 



/ „ , (21) 



and the lower dashed lines correspond to the positions of the chemical potentials ijS' 
and respectively, while the middle dashed line corresponds to the average value 
jj,. Parabolas (1) and (2) represent the dispersion law of electrons in wires 1 and 2 
respectively. The full circles correspond to the initial states of colliding electrons. 

Before the collision states la and 2a are occupied. The circle representing state la 
is below the dashed line, i.e. below the Fermi level /x. The circle 2a represents a state 
with p > 0, which is also occupied as the corresponding energy is below 

After the collision state lb is occupied. It is represented by a circle above the 
dashed line, which means that it has been free before the collision. In wire 2 state 2b 
with p < is also occupied. It is above i.e. it had been free before the transition. 

The width of the stripe between the two straight lines is eV. If the bottoms of the 
active subbands are well below the Fermi level the drag current should be proportional 
to the number of the occupied initial states as well as to the number of free final states. 

To calculate the drag current, one can recast the product of distribution functions 
in the collision term Eqs. (0) and (|8]) into the form 

V = 2sinh(eV72T)Q (22) 



where 



Q = exp exp - ,)f (e% ~ /x< + >) 

x /(effi - ai)/ - aiH) . (23) 
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For the drag current one gets 

. , feV\8e 5 mL ^ f°° „ f°° „ , gnn ,(k + k') _ 

' = — H^J^^io */„ <wS T+F ift (24) 

As above, one can conclude that the terms which give the main contribution to the 
drag current are those where \eiP (0) — e^/(0)| is smaller than or of the order of k B T or 
eV. We will assume that there is only one such difference (otherwise we would have got 
a sum of several terms of the same structure). 

As Q is a sharp function of k and k', one can take out of the integral all the slowly 
varying functions and get (the result is given for the general case where (jf- 1 * ^ // 2 )) 

eV e nn > eV e nn > 



7 ,I sinh I _fM **bT 2k B T 4k B T 2k B T 



Ak B T 2k B T J \4:k B T 2k B T 

where 

J ° = — Vnn '- (26) 

Herein' = e£\o) - e$ {0) . 

For eV <C k B T Eq. (p5|) turns into Eq. (|l6l). Let us consider the opposite case 
eV k B T. One gets for the drag current 



B 



eV 



This result is nonvanishing only if \e nn ,\ < eV/2. 

In this Section we have discussed a Fermi liquid theory of the Coulomb drag current 
in a quantum wire brought about by a current in a nearby parallel quantum wire. 
A ballistic transport in both quantum wires is assumed. The drag current Id as a 
function of the wire widths comprises one or several spikes; the position of each spike 
is determined by a coincidence of a pair of levels of transverse quantization, e n (0) and 
£ n '(0) in both wires. 



Figure 2. Schematic representation (for /Z 1 ) = /i^) of simultaneous transitions due 
to the interaction between electrons of the two wires for eV ^S> kBT. Circles o and • 
represent the initially unoccupied and occupied states respectively 



Figure 3. I /Id 1S plotted (for fj,^ = /i' 2 ' = fi) as a function of W\/W<i where 
the width of wire 1 is controlled through gate voltage (fi =14meV, T —1 Kelvin, 
W2 = 42 nm, L— 1 fixa, k = 13 and the spacing between wires is 50 ran). 
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3. Luttinger Liquid Theory of Coulomb Drag 

In ID systems e-e interaction gives rise to electronic correlations that are believed to 
destroy the Fermi liquid. Instead, a different state is generated that is usually described 
as a Luttinger liquid [[!7], [L9| (for reviews see e.g. ^(| |38], [3^, (|(J). It is therefore 



not surprising that in ID systems e-e interaction affects the drag in a different way 
than in two- or three dimensional systems. Indeed, in ID systems interaction strongly 
enhances the effect, getting stronger the lower the temperature. As a result, the positive 
temperature characteristic of the drag resistance from Fermi liquid theory can become 
a negative one. For sufficiently long wires the drag resistance becomes exponentially 
large at low temperatures. 



This chapter reviews in the main part the works [15, Id, 17|, and is organized 



as follows: Sec. pTT] introduces bosonic variables as the appropriate language for the 



discussion to follow. In Sec. |3.2| the renormalization group method is employed in order 
to show that and in which way the drag becomes enhanced by electron correlations. This 
consideration will also clarify some relevant energy and length scales of the problem. 



Sec. [O] elaborates on the influence of the electron spin on the drag, while Sec. [3]4| deals 
with non-linearities and asymmetric double-wires. The last section of this chapter, 
Sec. |3.5|, briefly discusses the drag in a system with a finite region of interaction. 



3.1. Bosonic variables 

For treating interactions it is convenient to describe excitations of the many-electron 
system in terms of collective coordinates: for example by the displacement ip(x) of 
electrons. They are normalized in such a way that density and current fluctuations are 
given by d x ip(x) = — ^ l /^l(n(x) — n ) and d t (p(x) = s/irl(x)/e. Rewriting the Hamilton 
of an interacting ID electron gas in <p(x) and its canonical conjugated field U(x) yields 
the Hamiltonian of an elastic string 

H=\j Mn 2 + i(M 2 . (28) 

The stiffness or interaction parameter K and the velocity v are determined by the 
parameters of the electronic system. For non- interacting electrons K = 1, v = vp, 
while for a system with repulsive interaction < K < 1 and v ~ Vp/K. The solutions 
(p(x, t) of the Hamiltonian (^) are ID waves with wave velocity v. In the limit of strong 
interactions K <C 1, these solutions correspond to the plasma oscillations of the electron 
density. Contrary to the underlying fermionic operators, the fields ip and n obey bosonic 
commutation relations. The substitution of the former by the latter is therefore known 



under the name of "bosonization" [19, 2C, Bq, pM 40 



Excitations of a double wire can be similarly described by the respective 
displacement fields <fi(x) and <p2{x) of each wire. Assuming a symmetrical system, the 
two eigenmodes of the density oscillations (at a given wavenumber q) are the symmetric 
mode (+), where the density in both wires oscillate in phase, and the anti-symmetric 
mode (— ), where the phases of the density oscillations differ by tt. A Transformation 
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to the corresponding displacement fields 0± = (ipi ±<^ 2 )/\/2 decouples the Hamiltonian 
into a symmetric and anti-symmetric part, 

H =H 1 + H 2 , 

H ± = V fjdx K ± U 2 + ^-(d x <P ± ) 2 . 
Each part has its own set of parameters. In good approximation 

K± = ( 1 + }i±3i_v^y"\ (29) 

\ TIVf TTVp J 

and v± = vf/K±, where Vq, Vq are the Fourier transforms of intra- wire and inter- 
wire interaction V(x) and V(x) for small momentum q — > 0. Vzk F is the intrawire 
backscattering strength (5q = 2kp) HI?]] . For the range of applicability of expression 



(p9|) see Refs. jy], |4^1 , where more precise estimates of the Luttinger-liquid parameters 
are given. 

So far interwire backscattering of electrons, where a large momentum of order 
5q = 2k p is exchanged, has not been taken into account. As pointed out in the previous 
chapter, this coupling however is essential for the drag and must be incorporated in 
the description. Fortunately, this can be also done in terms of the displacement field, 
leading to 

E 2 r 

H b = X^- / tfecos(V87T0-) ■ (30) 

TTVp J 

The energy Eq is of order of the Fermi energy, and the dimensionless coupling A is given 
by 

A = -^, (31) 

Note that the symmetric and anti-symmetric modes remain still decoupled. There is no 
corresponding term in the intra-wire interaction. The reason is that the backscattering 
within a wire appears as the exchange part of the forward scattering (5q — * 0), and 
therefore can be absorbed in the parameters K± and v ± (cf. Eq. ( ^9|) ) . 

The origin of the backscattering Hamiltonian Hb becomes clear in the limit of 
strong repulsive intra-wire interaction. In this case the electrons of each wire form 
well-correlated states with charge densities periodic in 2ir/kF- Accordingly, their local 
interaction energy is 27r-periodic in the relative displacement (si — S2)kp = V8n(j)- The 
integral over cos v^87T0- (x) with an appropriate prefactor therefore gives to first order 
the corresponding part of the total energy. 



3.2. Drag 

The backscattering Hamiltonian (|30|) is of sine-Gordon type, and allows for an intuitive 
understanding of the drag in the case of large couplings A. Suppose that the total energy 
is dominated by Ht>, the system minimizes its energy by fixing the field </>_ to a value 
\^8n(po = 7r + 2Tcm, m an integer number. Accordingly, the relative displacement of 



Coulomb drag between ballistic one- dimensional electron systems 



12 



electrons in wire 1 and 2 is constant in time and space. This means that two interlocked 
charge density waves have been formed, such that a current in one wire is necessarily 
accompanied by an equally large current in the second wire. In this ideal situation the 
drag is absolute |L5] . 



What happens to the drag if the situation is not that ideal is the issue of this 
section. Considerable insight with a minimum of calculation effort will be gained by 
making use of the concept of renormalization. 

Neglecting for a while all interactions except for the inter-wire backscattering, the 
double wire system can be viewed as a pair of uncorrelated ID Fermi liquids. As 
explained in the previous section, the inter-wire backscattering coupling then causes a 
drag resistivity pn = Rd / L, where L is the length of the drag wire, linear in temperature 
and proportional to A 2 , 

p D « po\ 2 T/E (32) 

(cf. Eq. [16] in the limit of T ^> Ae n (0)). To first order the drag only depends on the 
direct backward scattering part of the interaction. However, higher order contributions 
to pz> include inter- and also intra- wire forward scattering. In ID these higher order 
contributions are crucial and must be taken into account. The renormalization group 
theory does the job quite elegantly by successively integrating out high energy degrees 
of freedom down to an energy scale E < E . As a result, the original ("bare") couplings 
A become renormalized to ^-dependent couplings K_(E) and X(E). The energy 
scale at which the renormalization procedure has to be stopped can be given either 
by the temperature, the system size, or even by the coupling X(E) itself, depending 
on the circumstances. The net effect of higher order processes on the drag between 
a pair of weakly coupled wire can be summarized by replacing - for example in Eq. 
([32]) - the bare coupling A by a renormalized energy dependent ("running") coupling 
constant X(E). Further, good approximations for relevant energy scales can be easily 
extracted from the renormalization procedure. This should be enough motivation for a 
short excursion to the renormalization flow of the sine-Gordon model. 



3.2.1. Renormalization flow The flow is well-known from a closely related problem, 
that of an interacting spin-1/2 electron liquid ||44| . For small couplings A 1 it is 
described by the differential equations 

£ = (2-2** ^h-^* 2 , (33, 

where t denotes the negative logarithm of the rescaled energy, t = In Eq/E. (The 
subscript "— " is suppressed henceforth.) Fig. [| shows schematically the flow in a 
K — A-diagram. Each point represents a system characterized by the parameters K 
and A. Under renormalization the system develops according to the stream lines in the 
parameter space. The arrows indicate the direction of decreasing energy scale. 

The main feature of the flow is the so-called Kosterlitz-Thouless transition: systems 
with parameter below the line K — 1 + A renormalize towards weaker backward 
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Figure 4. The RG-flow of a double wire system of spin-less electrons. Point 
A corresponds to the bare couplings of a double wire with a rather large inter-wire 
distance d^$> Xp, point B to wires with narrow spacing d< Af. Point C corresponds 
to a single spin-l/2-liquid. 

scattering A, systems above or left of this line renormalize towards larger A. Systems 
right on the transition line flow to the point at K = 1 and A = 0, which represents a 
non-interacting Fermi liquid. 

Whether a system develops to weaker or stronger couplings A obviously depends on 
the location of the bare coupling (subscript "0"), the initial points of the renormalization 
trajectories. For systems with symmetrical interaction, V(x) = V(x), as it is the case for 
a spin-1/2 electron liquid, one finds for small V^kp that K = 1 + Ao, i.e. the initial point 
lies exactly on the transition line. These systems (marginally) renormalize towards the 
non-interacting Fermi-point. This is indeed the expected behaviour for the spin- mode 
of an electron liquid. (The spin-mode corresponds to the anti-symmetric mode (— ).) If 
one could measure the relative drag of spin-up and spin-down electrons, one would find 
that due to the additional temperature dependence in A(T), the drag resistivity decays 
with temperature even faster than the naively expected linear behaviour. 

Interestingly, the situation is completely different for a real double-wire system. 
Due to the spatial separation of the two wires, the intra- wire interaction V(x) always 
exceeds the inter- wire interaction V(x). Inspection of Eq.s (p9| ) and (|3~lD reveals that 
as a consequence the initial point is always above or left of the transition line, where A 
renormalizes to higher values. Therefore here the backscattering coupling A, although 
usually much smaller than in the previous case, increases with decreasing energy scale 
or temperature. Hence, for low temperatures (and long wires, see below) the drag 
resistivity will be always larger than predicted by Fermi theory. This will be made more 
quantitatively in the next paragraph. 

3.2.2. Temperature dependence of the drag Before rushing into a discussion of the 
various kinds of regimes with different types of temperature dependencies, it is advisable 
to clarify the relevant energy and length scales. There is E ~ Ep, the largest energy 
scale, the temperature Tl, at which the thermal wavelength becomes of the order of 
the system size, and the actual temperature T. The corresponding length scales are 
the Fermi wavelength A^ ~ Eq/v, the system size L, and the thermal wavelength 
Ljn = v/T. Less obviously, a fourth energy and length scale are given by the sine- 
Gordon Hamiltonian if_ + Hf,: the energy (mass) M of a soliton and its width L s . The 
soliton mass M coincides with the energy scale at which the renormalization procedure 
breaks down, the soliton width L s is the corresponding length scale. The relative order 
of these scales classify several different regimes. 

For high temperatures T > Tl, M, the renormalization of A is terminated by T. 
Thus, A = A(T). In this case the temperature dependence of A(T) can be approximately 
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determined by integrating the flow equations (|33|), leading to 

(rp \ 2K— 2 
Y Q ) ■ ( 34 ) 

This result inserted in Eq. (|32|) gives the temperature dependence, 

Pd~PoAo(— J , (35) 

valid for Tl,M <C T <C £q. The e-e interaction changes via the renormalization of 
the backscattering coupling, the linear temperature dependence known from the ID 
Fermi-liquid to a temperature scaling with an interaction dependent power x — 4if — 3. 
For sufficiently strong interaction K can assume values below 3/4. Then the power \ 
becomes negative and the drag increases with decreasing temperature. For vanishing 
interaction, K = 1, Eq. (|3~5"D goes over to the linear behaviour of the Fermi-liquid. 

Lowering the temperature below M or Tl, two scenarios are possible: if the wire 
is sufficiently long, L 3> W, at a temperature T ~ M the system eventually enters the 
strongly coupled regime, where A(T) ~ 1. For short wires, L ~ L s , this regime can 
not be reached. Here the renormalization halts at a temperature T ~ Tl, where the 
thermal wavelength is of order the system size. In this case even at low temperatures 
the systems are weakly coupled, \(T L ) <C 1. 

In the strongly coupled regime the energy is dominated by the backscattering term 
Hb, giving rise to an almost absolute drag. Deviations from this ideal drag correspond 
to processes where the relative displacement \fSTx<p_ slips from one global minimum 
position, say at 7r, to a neighbouring one at — it or 37r. At temperatures T < M these 
processes are enabled by thermally activated solitons moving along the wire. As a result 
the drag resistivity shows for T < M an activated behaviour, 

Pd(T) ~ p e M/T • 

This behaviour changes again when the temperature falls below Tl (for the strong 
coupling regime considered, Tl < M). It has been shown [irBJ that then the drag 
decreases linearly with temperature, due to the set-in of coherent soliton-tunneling. At 
even lower temperatures T < ^/T L M exp —M/T L the drag decreases with T 2 |I"6| . 

In the weakly coupled regime, the drag resistance decays oc T 2 as the temperature 
drops below Tl- This can be understood as follows: having renormalized down to an 
energy Tl, the original Fermi wavelength length Xp ~ v/E has become enlarged to a 
rescaled wavelength Xf{Tl) ~ v/Tl = L. Hence, the wires of length L have effectively 
shrunk down to pointlike constrictions connecting electronic reservoirs on either side. 
The drag in this situation is equivalent to the drag in a pair of ID Fermi-liquids (K = 1) 
that interact over a short length ~ Lp only. As it will become clear in the Sec. |3.5| , the 
drag is then proportional to T 2 (cf. Sec. |3~5| , K — 1). 

The temperature scale T* at which the system enters the strongly coupled regime 
is given by the soliton mass T* = M. Estimating it by A(T*) ~ 1 with the approximate 
expression (|34|) yields 

T* ~ E \^ . (36) 
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3.3. Electron Spin 

For comparison to experiments the treatment of the electron spin is mandatory. To this 
end one can introduce bosonic fields ip c / s that are related to the charge/spin density 
n c = nj ± ni in the same way as above ip is related to the density n of spinless 
particles. For a double wire, this results in a total of four modes: symmetric and 
anti-symmetric charge modes (c+ and c— ) as before, and, additionally, symmetric and 
anti-symmetric spin modes (s+ and s— ). Each mode is again described by a quadratic 
Hamiltonian of the type (|28|) with corresponding interaction parameters K C ±,K S ±, etc. 
. The neutral spin-modes are not affected by the interaction, wherefore K s ± = 1 
and v s ± = vf- Nevertheless, despite their neutrality, the spin- modes weakly couple 
to the anti-symmetric charge mode c— via backscattering processes. Hence, the drag is 
influenced by the spin degree of freedom []T7| . 



A quantitative analysis of the weakly coupled regime can again be done by making 
use of renormalization along the lines described above. The main results are summarized 
below: In the presence of spin the interwire backscattering coupling scales towards 
stronger couplings. However, fluctuations in the neutral spin-modes moderate the 
enhancement due to the interactions. This is reflected by an effective interaction 
parameter 

K c ^ + K s± K c ^ + 1 

Keff = — 2 — = 

which is closer to the non- interacting value 1 than the original K c _ . The parameter K c _ 
is given by 

-1/2 

Kc _ ~ ( i + 2 ^ ^ _ -_4kf_ i (37) 



Vo-Vo V^, 

TTVf TTVf 

As a result, in the weakly coupled regime the drag resistance scales with temperature 
as 

PD « PoAq ( -gr 

(cf. Eq. (|35D). The cross-over temperature T* turns out to be approximately 

Comparison with Eq. ( j36|) reveals again the moderating effect of the spin. If two systems 
have similar interaction constants K_ m K c _, but one is spin-polarized while the other 
is not, their respective cross-over temperatures and lengths are related by 

T*\ 2 ( T*\ ( Xf\ 2 ( A f 



0/ pol. J un—pol. V ^* J pol. V / un—pol. 

Since T*/Eq ~ \f/L* is usually a small number, the cross-over temperature of the spin- 
unpolarized system is by orders of magnitudes smaller than the one of a comparable 
spin-polarized system. 
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3.4- Non-linear drag and mismatching Fermi-momenta 

So far our considerations were confined to the linear regime (/ — > 0) of a symmetrical 
double wire system. This section extends the discussion to both the non-linear regime, 
and to systems with a misfit in the Fermi momenta, 5k = kpi — kp 2 ^ 0. 

It is again useful to look first at the associated energies. A finite current I in the 
active wire defines an energy Q = I/e, and the energy associated to the misfit is of 
course A = vSk. Non-linearities of the drag voltage Vd in the current /, or an effect of 
the misfit 8k will be significant only if the corresponding energies \Q\ or |A| exceed T, 
T L and M. 

In the weakly coupled regime both cases can be analyzed by perturbative methods. 
Finite currents and a non-vanishing 5k can be treated by a transformation <p(x, t) — > 
4>(x,t) + Ax/v + Qt. The term linear in x describes the density difference (d x <f> oc 
(jii — 712)), the term linear in t corresponds to a Galilei boost of the active wire. 
Accordingly, Hb becomes 

E 2 f 

H b = A — — / dx cos V8tt(4> + Ax/v + Qt) . 

7TV F J 

It is possible to derive a closed expression for the drag voltage Vd as a response to this 



perturbation |13 |. It is valid for arbitrary ratios f2/A, and can be written as 



2\2 / rp \ 4K 



eV D „E%K ( T 



C=^M — {A(Q - A)B(£l + A) 
v \E J 

+ A(Q + A)B({1 - A)} . (38) 

C is a numerical constant of order unity, A and B are temperature dependent functions, 
given by 



/E ( 1 T s "\ 

ds i sin(— s) I 7r( — h i) sinh — — I 
E \ s ttE J 



-2K 



and a similar expression with cos instead of isin for B. This expression holds also for 
vanishing A, Q, where it leads to the result (|35D . For current and 5k large compared to 
temperature, Q, A 3> T, Eq. (p8|) reduces to |15| 



^ = const. \ 2 (n 2 - A 2 ) 2K -\ for |A| < |0| 

= otherwise . 

For non- vanishing A the drag vanishes as long as the current is below a threshold value. 
For larger currents the voltage shows a powerlaw dependence on the current. A thorough 
discussion of the non-linear I — V characteristic can be found in [JT5[| . 

Actually, at finite temperatures the drag does not vanish completely for |A| > |0|, 
rather it shows an activated behaviour as like in the case of ID Fermi liquids. This can 
be made more explicit jEJ. Evaluating expression ( P5| ) in the limit Q —>■ at finite A 
results in a drag resistivity 

Pa,t = PoA 2 ( — ) F 2K (A/T), (39) 
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where F 2K is an interaction dependent function defined by 

F 2K {e) = N 2K (e)^(-s) + N 2K (-s)^(e) , 

N 2K {e) = lim J ds e is£/7r (d+i) sinh s) 

F 2 k is a continuous function. It decays exponentially at large positive arguments, 
F 2 k{z) ~ exp — e, and behaves algebraically for large negative arguments, F 2K (e) ~ 
|e| . The result obtained for non-interacting ID Fermi-liquids is recovered by putting 
K = 1. In this case 

ir e £ — 1 

such that the expression (|39| ) corresponds to Eq. ([U]). While N 2 (e) bears some 
resemblance to the Bose distribution, for the special interaction parameter K = 1/2 
one obtains exactly the Fermi function 

1 



NAs) 



e £ + l 

For general parameter K an analytical expression for N 2 k is lacking. It is an open 
question whether the functions N 2 k are related to the exclusion statistics of fractional 
excitations in the Luttinger liquids. 



3.5. Finite interaction region 

Double wires that interact only over a region of finite length I < L have been also 
investigated |12|, [13|]. For temperatures T < v/l this problem can be mapped to 
the classical problem of a Luttinger liquid with a single impurity. Qualitatively, 
these systems behave similar to those considered in the previous sections. In the 
weakly coupled regime the drag scales with temperature with an interaction dependent 
exponent, 4K — 2. A strongly coupled regime with almost absolute drag at zero 
temperature exists as well. However, it is reached only for sufficiently strong interaction 
K < 1/2. For K > 1/2 the interwire backscattering coupling renormalizes to weaker 
couplings, such that the drag vanishes for T — >• 0. 



4. Experimental search for ID Coulomb drag 

Although a fair amount of theoretical work has been available on Coulomb drag between 
ID electron systems, there has been a conspicuous absence of experimental work. This 
may be attributed to two difficulties encountered in measuring the ID drag. First, since 
it is a very small effect, the drag voltage usually has a very small magnitude and must 
be clearly distinguished from spurious signals. Second, and perhaps the major difficulty, 
has been the difficulty in creating parallel, electrically isolated, quantum wires with a 
spatial separation large enough to completely suppress interwire, while small enough to 
give a drag voltage of a reasonable magnitude. It was only recently that Debray et al |J 
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reported the first experimental observation of Coulomb drag between ballistic quantum 
wires. The same authors later published a more comprehensive experimental work |J 
on the subject. Work along the same lines has lately been reported by Yamamoto et 
al [0. In the following, we give a brief outline of the reported experimental work and 
an analytical discussion of the results in the framework of the Fermi and the Luttinger 
liquid theory as discussed in Sees. II and III. 

4-1. A. Experimental techniques for dual-wire sample realization 

The samples used for ID Coulomb drag measurements consisted of two electrically 
isolated, parallel quantum wires, with a small spatial separation. Such samples were 
fabricated from AlGaAs/GaAs heterostructures with a high-mobility (=10 6 cm 2 /Vs) 
two-dimensional electron gas (2DEG) at the interface. The dual-wire samples 
were fabricated using high-resolution electron beam lithography, combined with deep 
chemical etching. The samples used so far were made in a planar geometry by depletion 
of a single 2DEG layer by three surface Schottky gates || |5], [7J deposited on the 
heterostructure wafer. Figure |5] gives a schematic top view of the planar device and 
the scanning electron micrograph of a typical device ||. U, M, and L are surface 
Schottky gates. 



Figure 5. (a) Schematic top view of a planar Coulomb drag device. U, M, and L are 
surface Schottky gates, (b) A scanning electron micrograph of a typical device with 
middle gate width of 50 nm. 

Dual-wire samples for drag measurements can also be made in a vertical geometry 
H from two vertically stacked quantum well (QW) structures with a 2DEG in each well 
(such samples have not yet been used). The advantage of the planar geometry is that the 
interwire separation can be changed in-situ by changing the bias voltage of the central 
gate M. The disadvantage is that the narrow central gate creates a soft lateral potential 
barrier and to prevent tunneling between the wires the width of this barrier has to be 
of a sufficient magnitude, which sets a limit to the minimum interwire distance that can 
be used without electron tunneling interfering. Samples with a vertical geometry have 
been widely used for studying Coulomb drag between 2D electron layers [2]. The main 
advantage of the vertical geometry is that very small interwire separation (the barrier 
width) can be obtained without tunneling between wires. Since the magnitude of the 
drag is expected to decrease exponentially with interwire separation, one can expect 
to observe enhanced drag with the vertical samples because of the smaller separation 
that can be achieved with such samples. The major disadvantage is that the interwire 
separation cannot be changed in-situ, in contrast to the planar case. Also, it is not 
obvious that the widths of the two wires can be independently changed through the use 
of the mutually aligned top and bottom split gates. 
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4-2. B. Experimental observation of Coulomb drag 

In their work, Debray et al pi used a planar geometry of quantum wires, as shown in 
Fig. 5, of lithographic length L= 2 /mi with a middle gate M of lithographic width 50 
nm. The drag voltage Vd was measured with a drive voltage Vds i n the linear regime 
of ballistic electron transport as a function of the width of the drive wire by adjusting 
the bias voltage Vu, while the width of the drag wire was adjusted to have the Fermi 
level Ep just above the bottom of its lowest ID subband. An appropriate negative bias 
voltage Vm was applied to the middle gate to ensure total absence of interwire tunneling. 
Measurements were done in the absence of any such tunneling. In Fig. |6]a is shown the 
measured drag voltage Vd as a function of the width of the drive wire. The drag voltage 
is found to show peaks, which occur in the rising parts between the plateaus of the drive 
wire conductance. This suggests that they occur when the ID subbands of the wires are 
aligned and the Fermi wave vector kp is small. Measurements carried out in a magnetic 
field B = 0.86 Tesla perpendicular to the plane of the device shown in Fig.6(b) indicate 
identical behavior except that the magnitude of Vd is enhanced almost by a factor of 
three. 



Figure 6. The drag voltage Vd and drive current / as function of drive wire width 
at a drive voltage Vds = 300 fjY M . (a) In zero magnetic field with the upper wire as 
the drive wire, (b) In a magnetic field of 0.86 Tesla with the bottom wire as the drive 
wire. 

In order to have unambiguous evidence that the observed drag voltage Vd is indeed 
due to the Coulomb drag effect, the authors measured the dependence of Vd and Rd on 
the interwire separation and the temperature. Figure [7| shows the dependence of Vdm, 
the height of the first Vd peak of Fig. ||b and the corresponding Rd as function of the 
middle gate bias voltage Vm- The two quantum wires (Fig. 0) were spatially separated 
by an effective distance d due to the depletion by Vm of the 2DEG under the middle 
gate M. In the voltage range of interest, d was experimentally found to vary almost 
linearly with Vm according to 

d = d + a {V - V M ) , (40) 

where Vo is the value of Vm for which the 2DEG under M is just depleted and a gives 
the total spatial displacement of the two depletion edges of M with respect to its bias 
voltage, do is a constant for the same device and is nominally equal to the lithographic 
width of the gate M. One can change d by varying Vm- In Eq.(^), V and a were 
determined experimentally. The dependence of Rd on Vm was found to be exponential 
and can be described well by the relation, Rd oc e^ Vu , where (3 = 14.2(9)V _1 . 

The temperature dependence of Coulomb drag is a crucial feature that can be used 
to probe which one of the two theoretical models, the FL or the LL theory, constitutes a 
more appropriate description of ID Coulomb drag. Measurements, such as shown in Fig. 
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Figure 7. Dependence of the drag response on interwire separation d via the middle 
gate voltage Vm |§- (a) The maximum Vdm of the first drag peak of Fig. 6(b) as a 
function of Vm- (b) The natural logarithm of the corresponding drag resistance Rd 
as a function of Vm- The dotted line is a linear fit to the data points. 

H carried out in the temperature range 60 mK - 1.2K are shown in Fig. A decrease of 
Vd with increasing temperature was observed. The dependence on temperature of the 
drag resistance Rd corresponding to Vdm is shown in Fig. [| for both in the absence and 
presence of an applied magnetic field B. The temperature dependence can be described 
well by the power law, Rd oc T x , with x = -0.77(2) and -0.73(6) for B = and B = 0.86 
Tesla, respectively. It is interesting to note that the data points at temperatures lower 
than 180 mK, for zero field, and 300 mK, for nonzero field, fall below the power-law 
curve, indicating a suppression of the drag effect. 

Figure 8. The dependence of drag voltage on temperature ||. (a) The drag voltage 
Vd as a function of the width of the upper (drive) wire in zero magnetic field with 
Vds = 300 /iV at 70, 180, 300, 450, and 900 mK, corresponding to curves in order of 
decreasing peak height, (b) The same as in (a) but in a magnetic field of 0.86 tesla 
with V DS = 50 ^V at 60, 180, 300, 450, 900 mK, and 1.2K. 



Figure 9. The temperature dependence of drag resistance Rd corresponding to VnMof 
the first drag peak of Fig. 8 in zero field (a) and in a field of 0.86 Tesla (b) Q. Note 
that the data points at the low end of the temperature range fall below the power-law 
curve. 

Lately, using a lateral sample geometry , very similar to that shown in Fig. [5]a, 
Yamamoto et al |7| has reported the observation of Coulomb drag and the influence of 
an applied magnetic field on it. Their results corroborate those of Debray et al |5j ||. In 
their work, Yamamoto et al also reported the observation of a negative drag. However, 
since the negative drag was observed only when the drive wire was completely pinched 
off, it is highly questionable if the effect observed is due to Coulomb drag. 

4-3. C. Discussion 

The origin of the observed peaks in the drag voltage Vd (Fig. 6) can be understood when 
one considers Eqs. (16)-(21) of Sec. II. Since Vd is directly proportional to the drag 
current Id, Vd will show maxima whenever any two ID subband bottoms of the two 
wires line up and the Fermi wave vectors in the two wires are equal and small. As seen 
from Fig. 6, the occurrences of the drag peaks correspond to these conditions. The 
first peak in Vd occurs when the Fermi level is just above the bottoms of the lowest ID 
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subbands of both wires. Similarly, the second peak occurs when the lowest subband of 
the drag wire lines up with the second subband of the drive wire. Both the increase 
and the narrowing of the first drag peak in a magnetic field of 0.86 Tesla (Fig 6(b)) 
can be attributed to an increase of the density of states inlD the subbands due to the 
magnetic-field-induced enhancement of the electron effective mass. The reduction in the 
magnitude of the drag peak as we move away from the first peak toward higher values 
of Vu can be attributed to an increase in the effective interwire separation of the wires. 
This dependence is explained in detail later. Since we are mainly concerned with ID 
transport in the fundamental mode, we restrict our discussion here to the region of the 
first drag peak. 

To understand the dependence of drag on the interwire distance shown in Fig. 7, we 
note that the matrix element of the backscattering probability depends on the interwire 
distance d via the modified Bessel function KQ{2kpd) (Sec. II, Eq.(ll)), which is an 
exponential function of its argument for 2k pd » 1. The same dependence also results 
from the LL theory [29]. Under this condition, an exponential decrease of Rd with d is 
expected according to Rd oc exp (— Ak F d). This is consistent with the results of Fig. 7. 
Using the experimentally determined values (3 = 14.2 V -1 and Vo = -0.4V, a = 580 
nm V -1 , we find ^ = 6.1 x 10 6 m _1 . Surprisingly, this corresponds to a low density of 
about 8 electrons per 2/xm wire segment and a mean electron distance f « 250nm in the 
wire. When Vm is in the range from -0.7 to -0.8 V, we have (Eq. (fK]) d = 0.2 /j,m. This 
gives 2kpd = 3, so the approximation of 2fcp(i > 1 is reasonable. This exponential 
decrease of Rd with d also explains why the height of the drag voltage peaks in Fig. 6 
decreases so rapidly as Vu increases. An increase in Vu increases the width of the drive 
wire and hence d. The decrease of Rd for Vm > — 0.7V occurs due to tunneling of a 
considerable fraction of the current from the drive wire to the drag one, reducing the 
measured Rd- 

The experimental observed features of the drag effect discussed above, namely, 
the origin of the drag voltage peaks, the effect the magnetic field, and the interwire 
separation dependence, can all be understood in the framework of both the FL and the 
LL theory. It is the temperature dependence of the drag that is the crucial feature - it 
can be used to determine which one of the two theoretical models constitutes a more 
appropriate description of ID Coulomb drag observed under the given experimental 
conditions. The observed temperature dependence of Rd, shown in Fig.9, is in sharp 
contrast with the linear temperature dependence predicted by the FL theory (Eq. (19)). 
The unusual temperature dependence can not be attributed to a temperature induced 
modification of the wire conductance, since the latter is found to be almost unchanged 
over the temperature range of the measurements. A reduction of the interwire Coulomb 
coupling due to enhanced screening by the reservoirs and gates is very unlikely at such 
small temperatures. On the other hand, it's conceivable that a correlated LL behavior is 
established in the wires. Indeed, it is hardly surprising that the temperature dependence 
of Rd does not fit into a FL scenario, because for the experimental condition of the first 
drag peak the ratio r s of f and the Bohr radius as, r s — f/as ^ 26 is large. 
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The smallness of the drag resistance (Rd < 1000) in zero magnetic field indicates 
a weak interwire back scattering coupling. In this case, according to the LL model Rd 
should obey a power law as long as the thermal length Lt is well in between the wire 
length L = 2 and the mean electron distance f rs 250nm in the wire. For spin- 
unpolarized electrons, valid for the data shown in Fig. 9, the LL description predicts 
a power-law temperature dependence of Rd with exponent x = 2K C _ — 1 (Eq.(38)). 
The data shown in Fig. 9 indeed shows a power-law dependence of Rd on temperature 
with K c \ = 0.12. Let us see if the condition f < Lp < L is fulfilled in the experiment. 
Given a Fermi wavevector of kp ~ 6/zm _1 we find that Lp = hvp / K c _kpT is equal to 
the wire length L = 2 /im at a temperature Tp = 250 mK, and that Lt approaches 
f « 250nm at a temperature of about 2K. Here Vp/K c _ is the group velocity of the 
relative electron-density fluctuations and h/kpT is the quantum lifetime associated with 
the thermal energy ksT. This means that there is a narrow temperature range window in 
which a power-law temperature dependence of Rd might be expected, and it is observed 
experimentally. At temperatures below Tp, when L < L-r,the electron coming from the 
lead to the wire does not have time to accommodate itself to the LL liquid. This should 
result in a drag weaker than the power-law dependence. The experimental data of 
Fig. 9 are consistent with this analysis. At lower temperatures we do indeed observe a 
tendency to a weakening of the drag with respect to the power-law dependence. 

The negative power-law temperature dependence is not the only experimental 
feature that can not be understood in terms of the FL theory of Coulomb drag. The 
experimental value of Rd (Fig. 9a) at T = 60 mK is more than an order of magnitude 
larger than that given by FL theory (Eq. (19)). That the measured drag is larger could 
be explained by the interaction-normalized interwire backscattering probability, which 
should be larger than the bare one (Eq. (34)). 

Comparison of Figs. 9(a) and 9(b) show that the influence of the magnetic field 
on the temperature dependence is not significant. This may signify that Zeeman spin 
splitting at B < 1 Tesla is not important yet, otherwise the exponent x should change 
(Eq. (35)). Indeed, a clear signature of spin splitting in the measured conductance 
staircase was not observed at this field. The magnetic field, however, increases Vp for 
the same position of the Fermi level. This makes Lt larger at the same temperature 
compared to that in zero field. This can explain why in a magnetic field a deviation 
from the power-law dependence occurs at a higher temperature (Fig. 9(b)). 

We have interpreted above the experimental data in terms of Coulomb drag only. 
Considering the large interwire separation for which the drag measurements were made, 
one can not rule out the possibility of an acoustic phonon-mediated drag (PMD) 
contribution to the measured drag resistance. Recent theoretical work |45j j|6l on ID 
PMD based on Fermi liquid description predicts that PMD is negligible compared to 
Coulomb drag for 2kpd < 5. Also, for a dual-wire sample shown in Fig. 5, Rd should 
increase exponentially with temperature in the range 100 - 600 mK and does not decrease 
exponentially with interwire separation d. The data shown in Figs. 7 and 9 qualitatively 
contradict these predictions. This allows us to conclude that the PMD contribution, if 
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present at all, is insignificant. 



5. V. Future prospects 

It is quite obvious from the content of Sec. IV that substantial experimental work 
remains to be done to gain a comprehensive understanding of the physics of Coulomb 
drag between interacting ID electron systems and to explore the conditions under which 
such systems behave as a Fermi liquid or a Luttinger liquid. Since the measurement 
of the Coulomb drag also provides a new experimental tool to probe the LL state that 
can't be done from the measurement of the conductance alone, extensive experimental 
work on the subject is needed to put the LL model of interacting ID systems on a firm 
footing. Though the theory of Coulomb drag has considerably outstripped experimental 
work, many open questions need to be addressed in the theoretical well. 

On the experimental side, work should be focused on measurements that can 
distinguish between a LL and a FL state and can provide information about the existence 
and the nature of the LL state. This is an extremely important area for condensed matter 
physics. The few papers published so far, claiming to have observed a Luttinger liquid, 
have not been convincing. In this respect, it would be highly interesting to study the 
drag between spin polarized systems, since the LL theory predicts different exponents 
for spin polarized and unpolarized cases and manifestation of the spin effect should be 
quite different in the Fermi liquid and the Luttinger liquid state. Another interesting 
experimental possibility is to study drag when the wire length L falls below the thermal 
length Lt to investigate if the drag resistance Rd decays oc T 2 as predicted by LL 
liquid theory. When the number of electrons in the wires is very small (Sec. IV), one 
should expect relatively large fluctuations of the drag current or voltage, such as shot 
|47||[48l|, and possible reversal of the sign of drag leading to negative drag ||49|1 . 



noise 



Observation of this noise can also provide valuable information on correlated electron 



state. One could also envision a search for ID spin Coulomb drag [50|. Finally, it is also 



important to study acoustic phonon-mediated drag (PMD) ||5| [f|(J since under certain 
conditions it can be comparable to and even larger than the Coulomb drag. If such a 
PMD is present in the experimental measurements, one has to find ways to separate it 
from the Coulomb drag. 

The theory of Coulomb drag based on the LL model is far from mature and many 
open questions need to be addressed such as the effect of disorder, the influence of 
tunneling between the wires, etc.. Alhough the power-law temperature dependence of 
the drag resistance is a signature of the Luttinger liquid state, a careful analysis of 
various limiting cases based on the Fermi liquid approach should be carried out to make 
sure that under no circumstances it can give a similar temperature dependence. It 
is equally important to investigate the physical situations and interactions (within the 
wires and with the reservoirs) that favor transition of the Fermi liquid into the Luttinger 
liquid and vice versa. 
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